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ABSTRACT 

In hierarchical structure formation, merging of galaxies is frequent and known to dramatically affect their 
properties. To comprehend these interactions high-resolution simulations are indispensable because of the non- 
linear coupling between pc and Mpc scales. To this end, we present the first adaptive mesh refinement (AMR) 
simulation of two merging, low mass, initially gas-rich galaxies (1.8 x lO^^M© each), including star formation 
and feedback. With galaxies resolved by ~2 x 10^ total computational elements, we achieve unprecedented 
resolution of the multiphase interstellar medium, finding a widespread starburst in the merging galaxies via 
shock-induced star formation. The high dynamic range of AMR also allows us to follow the interplay between 
the galaxies and their embedding medium depicting how galactic outflows and a hot metal-rich halo form. 
These results demonstrate that AMR provides a powerful tool in understanding interacting galaxies. 
Subject headings: galaxies: interactions — galaxies: formation — galaxies: starburst — stars: formation 
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1. INTRODUCTION 

Decades of work have been devoted to the study of interact- 
ing and merging galaxies, as they play ess ential roles not only 
in shaping present-day galaxies (IToomreiri977i "merger hy- 
pothesis"), but a lso in constructing lar ge scale structures from 
the bottom up (IWhite & ReesI 119781 "hierarchical structure 
formation"). Because of the nonlinear coupling between pc 
(star forming regions) and Mpc scales (the distance at which 
tidal interactions occur) accurate numerical studies are imper- 
ative to comprehend the evolution of interacting galaxies. Al- 
though the morphology of merger remnants has been well re- 
pr oduced by N -body simulations since the pioneering work 
by lToomre & Toomre (1972), and various physical character- 
istics and merger-driven starbursts have been successfully an- 
alyzed with smoothed particle hydrodynamics (SPH) simu- 
lations jB arnes & Hernquist 1996; Mihos & Hernquist 1996; 
ISDringelet al. 2005a,b; Di Matteo et al. 2005; ' Cox et al.] 
|2006a,b; Mayer et al. 2007; Li et al. 2007; Saitoh et al. 2008), 
a complete, self-consistent simulation of galaxy mergers has 
not yet been perfected. 

First, SPH simulations tend to have coarse resolution in 
an interstellar medium (ISM), leading to the over-mixture of 
different gas phases ( Agertz et al. 2007; Tasker et al. 2008). 
Therefore, straightforward SPH simulations might have com- 
plications in realizing a multiphase medium, in capturing 
shock-induced star formation, in converting thermal feedback 
to a kinetic motion, and thus in showing how feedback makes 
a difference in galactic evolution self-consistently, though dif- 
ferent formulations and sub-resolution m odels alleviated the 
problem s (Marri & Whi te 2003: Springel ^Hernquistll200l 
lBarnesl l2004: Scannapi eco et al.ll2006h . 

Second, since gaseous halos and an intergalactic medium 
(IGM) have not been sufficiently resolved in SPH simula- 
tions, it is not easy to investigate the interplay between a 
galactic disk and a diffuse embedding medium. For instance, 
ICox et al.] (|2006a|) emphasized that a galactic halo should be 
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included to accurately study the galactic wind and the enrich- 
ment process powered by feedback and mergers. Yet because 
of their Lagrangian nature and smoothing scheme, SPH simu- 
lations might not have sufficient resolution to follow the evo- 
lution of the diffuse medium and the galactic outflow. 

In light of these needs, an adaptive mesh refinement 
(AMR) technique potentially provides a uniquely useful 
tool to address these issues, allowing us to realize a self- 
consistent, high-resolution galaxy merger simulation. As 
proven by an increasing number of groups (Tassis et al. 20031 
2008; Kravtsov & Gnedin] 120051: ICeverino & KlvpinI \200% 
iDubois & Teyssier 2008) AMR simulations have been highly 
successful in resolving the detailed structure of galactic evo- 
lution. In order to make use of the advantages of AMR such 
as the high dynamic range and r eliable shock resolution, we 
utilize the AMR code enzo (Brva n et al. 1120011: lO'Shea et aP 
I2OO 4J). In this Letter, we focus on the the first of its kind AMR 
simulation of two merging, low mass, initially gas-rich galax- 
ies, including star formation and feedback, with emphases on 
shock-induced star formation and the hot gas outflows. 

2. METHODOLOGY AND NUMERICAL SIMULATIONS 

2.1. Data Conversion Pipeline 

We developed a data conversion method which c onverts a 
galac tic N-body dataset of GalactlCS (Kuiiken & Dubinsl; 
Il995h to an SPH dataset for Gadget (Springel eTaP 1200 Ih 
and then to an adaptive mesh for enzo employing Delaunay 
tessellation onto an oct-tree structure. Using the particle data 
of a galactic sized halo with both the dark matter and the gas 
[Kim et al. ( 2008) demonstrated the compatibility of the initial 
N-body dataset and the adaptive mesh produced through the 
pipeline. A suite of functionality checks finds very satisfac- 
tory results enabling us to study galaxy evolution with AMR. 

2.2. Simulation Code 

The high-resolution Eulerian AMR code enzo captures the 
gravitational collapse of t urbulent fra gmentation with very 
high spatial resolution (^e.g. lWise et al.| [2008) and attains mul- 
tiphase gas dynamics in the ISM as it sharply resolves shocks 
and phase b oundaries (Tassis etal. 2003; Slyzetal. 2001: 
lAgertz et"aD 120071: iTasker et al. 2008). Enzo also contains 
all relevant physics previously discussed in simul ating galaxy 
evolution processes (ITasker & Brvanll2006ll2008h . 
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Fig. 1 . — The global K-S relation for an isolated galaxy: time variation of 
the relationship between global SFR and gas surface density. Each data point 
represents a different epoch, equally spaced in 5 Gyr. The solid line is the 
best fit for simulated data of Zsfr > lO~^M0yr~^kpc~^, and the dashed line 
for ZsFR > IQ-^Moyr-^kpc-^. The dotted line is from lKennicutl <1998h . 

We employ the ZEUS hydrodynamics module included in 
enzo to evolve the gas. Ra diative cooling is used by adopt- 
ing Sarazin & White! (I l987|) to follow the equilibrium cooling 
function down to 10^ K, and lRosen & BregmanI (Il995|) further 
down to 300 K. This treatment will ensure a thin galactic disk 
forms by being cooled below 10^ K, the approximate Tyir of 
the ISM in a galactic disk. The cutoff at 300 K roughly models 
the temperature floor provided by nonthermal pressure such as 
cosmic rays and magnetic fields (iRosen & Bregmanlll995|) . 

Galaxies are placed in a box of 4 Mpc on a side to ensure 
enough space for galactic tidal interactions and to reduce any 
boundary effect. The top grid of 128^ cells is allowed to recur- 
sively refine up to 13 levels based on the baryonic mass and 
the dark matter mass in each cell, achieving 3.8 pc resolution 
in the ISM. This value is in accord with the Jeans length for 
a dense gas clump of ^ = 10^ cm~^, at which a correspond- 
ing Jeans mass of 2 x lO^M© collapses to form a star particle. 
In this way, merging galaxies are resolved with ~2 x 10^ to- 
tal computational elements, surpassing any numerical studies 
conducted thus far on galaxy mergers includi ng gas. 

Our star formation criteria are based on ICen & Qstriked 
([l992) with several important modifications. A cell of size Ax 
produces a star particle of m* = epgasAx^ (e=0.5, a star forma- 
tion efficiency) when (i) the gas density exceeds ^thres = 10^ 
cm~^, the flow is converging, (Hi) the cooling time is 
shorter than the dynamical time, and ( iv) the particle produced 
has at least lO^M^. We do not impose any stochastic star for- 
mation unlike Tasker & BryanI (|^06) or Stins on et al.l (|2006). 
With these revisions, our criteria guarantee that a star parti- 
cle forms before an unphysically large mass begins to accrete 
onto any unresolved dense gas clump. 

The energy loss by radiative cooling can be replenished 
by thermal stellar feedback. For each star particle, 5 x 10~^ 
of its rest mass energy and 25% of its mass are returned to 
the gas over the dynamical time of the particle. This cor- 
responds to 10^^ ergs per every llOM© deposited as stel- 
lar mass and represents various types of feedback such as 
proto stellar outflows (Li & Nakamura 200 6). photoionization 
(iMcK ee 1989), stellar winds, and Type II supernovae explo- 
sions (Tas ker &Brvanll2006b . This feedback heats -10%© 
in a <10 pc cell up to ^10^ K, but a multiphase medium is 
naturally established because the cooling time of these hot 
cells is always much longer than the sound crossing time. 

2.3. Initial Conditions 

The individual galaxy progenitor we modeled has a total 
mass of 1 .8 X lO^^M© with 10% in gas (Ryir = 65 kpc at z=0). 
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Fig. 2. — The global SFR as a function of time during the galaxy merger 
run and the isolated galaxy formation run (twice the value). The separation 
between the centers of two galaxies is also displayed. 

Because we generate gas grids by splitting 10^ collisionless 
particles in N-body data with the same density profile and ve- 
locity dispersion, the gas will virialize to the desired Tytr of 
the galactic halo automatically. A spin parameter A =0.05 5 
is given to cause the progenitor to form a disk galaxy with 
a gaseous halo within a few hundred Myrs. In addition, this 
galaxy progenitor is bathed in a warm (10^ K) diffuse back- 
ground IGM; an initial metallicity of lO^^Z© is also set up 
everywhere to follow the metallicity evolution. For a merger 
simulation, two identical galaxy progenitors are separated by 
100 kpc and set on a prograde hyperbolic (^=1.1) coplanar 
collision course with a pericentric distance of 4 kpc. The ini- 
tial separation is large enough to form individual galaxies be- 
fore the first passage, and to observe the collision interface 
between the two gaseous halos. 

3. RESULTS 

3.1. Properties of An Isolated Galaxy Model 

We first examine how well our isolated galaxy formation 
simulation fits the global Kennicutt- Schmidt (K-S) relation 
between global star formation rat e (SFR) and gas surface 
density, namely Lsfr ^gas^'^ ( Kennicutt 1998). To cal- 
culate both densities, we select a disk of radius 2.8 kpc 
so that 95% of created stars is contained in it after 5 Gyr; 
the gas surface density is averaged over the cells where the 
gas density exceeds n = 4.0 x 10~^cm~^. Figure [T] shows 
our star formation criteria and feedback correctly match 
the observed K-S relation, as other simulation works did 
(e.g. [Robertson & Kravtsovl [2008). The closest match to 
the K-S relation occurs when we restrict the fit to Lsfr > 
lO^'^Moyr^^kpc"^, which happens m ostly in the first 2 Gyr, 
following an observational cutoff as in lKennicuttI (I l998|) . 

3.2. Star Formation History in A Galaxy Merger 

The global SFR of the merger simulation is displayed in 
Figure [21 It presents the initial stellar disk formation for each 
galaxy in the first ~0.6 Gyr, and several merger-driven star- 
bursts afterwards, notably when two galaxies first encounter 
(~0.8 Gyr) and when they finally coalesce to form one galaxy 
(~4.7 Gyr). A low SFR between these two bursts confirms the 
regulated star formation by stellar feedback. 

Snapshots of the merger sequence at four different epochs 
are compiled in Figure O The top row shows the density- 
weighted projection of density, in which irregular gas fila- 
ments, bridges and rings are formed by the compression of 
gas and turbulence. The middle row depicts the temperature 
sliced at the collisional plane, where cold gas clumps and hot 
supernovae bubbles coexist side by side forming a complex, 
yet well-resolved multiphase medium. It also reveals how hot 
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Fig. 3. — Density-weighted projection of density (top), temperature sliced at the orbital plane (middle), an d stellar distributions color ed by creation time 
(bottom) in the central 40 kpc, after 0.8, 1.3, 4.5, and 5.7 Gyr. High-resolution images and movies are available at |http://www. slac. Stanford. edu/^mornkr/] 



supernovae bubbles propagate through the diffuse embedding 
medium of the halo and the IGM. In the bottom row of stellar 
distributions colored by creation time, both merger-induced 
nuclear starbursts and shock-induced widespread starbursts 
are noticeable. Because of the finer resolution in the ISM, 
it is easier to resolve local dense clumps driven by shocks and 
the ensuing star formation. In contrast, SPH simulations often 
report predominantly nuclear starbursts (Barnes 2004) . 

3.3. Gas Outflows and Formation of A Hot Gaseous Halo 

The evolution of the stellar and gas mass in the central 200 
kpc box is plotted in Figure (H The gas expulsion via stellar 
feedback and galactic interaction is pronounced as more than 
90% of the gas has been expelled in the first 4 Gyr. This gas 
eventually escapes the gravitational potential of the system or 
has not had enough time to fall back onto the galaxies. This 
massive gas depletion is prominent especially in low mass 
mergers because of the shallow gravitational well. As for the 
merger remnant at 5.7 Gyr, the remaining gas mass is ~35% 
of the stellar mass and still decreasing rapidly. The amount 
of cold gas (T < 10^ K) available for future star formation is 
only <1% of the stellar mass, depicting how star formation is 
quenched by feedback heating and gas expulsion. 

The gas disrupted by galactic interaction and heated by 
feedback creates a galactic wind of >200 km/s reaching as 
far as 1 Mpc from the simulated merging galaxies. This hot 
metal-loaded outflow is responsible for building the gaseous 
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Fig. 4. — The evolution of the stellar mass (solid line for the merger run 
and dotted line for twice the value of the isolated galaxy) and the gas mass 
remaining in the central 200 kpc cube (short dashed line). Cold (T < 10^ 
K) and hot (T > 10^ K) gas masses are also shown (long dashed line and 
dot-dash line, respectively). 



halo around galaxies as well as enriching some regions of the 
IGM up to a supersolar metallicity. As a result, a hot metal- 
rich halo is generated (p ^ 10~^^gcm~^, T ~ 10^-^ K) and 
sustained by contin uous stellar feedb ack, as suggested by an- 
alytic models (e.g. iTang et al.|[20Q8|) . Although the galactic 
outflows and the halo are very diffuse, their evolution is eas- 
ily followed in AMR, as can be clearly seen in Figure \5\ of 
the joint probability distribution functions (PDFs) on density- 
temperature planes. It also illustrates the wide range of densi- 
ties and temperatures that are followed here. 
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Fig. 5. — Joint density-temperature PDFs, colored by metallicity, for a 400 
kpc sphere centered on the galaxy, after 0.8 Gyr (left) and 5.7 Gyr (right). Star 
forming regions, supernovae bubbles, gas outflows, and the halo are pointed 
out, proving the wide range of densities and temperatures followed here. 

4. DISCUSSION AND SUMMARY 

Our simulation, for the first time, followed the self- 
consistent evolution of low mass merging galaxies with AMR 
at unprecedented resolution. Our findings are as follows. 

First, as AMR naturally establishes a multiphase medium 
without any sub-resolution model, we have captured shock- 
induced star formation that occ urs when merging galax- 
ies co mpress the intervening gas (|Barnesll2QQ4l: [Saitoh et al." 
|2Q08|) . The well-resolved shocks trigger a widespread 
starbu rst, in accord with observations (e.g. ISchweizej 
2QQ6h. Further, the ove rcooling problem is absent as in 



Ceverino & KlypinI (l2QQ7l) because the multiphase medium is 



resolved by <10 pc cells, and the thermal feedback is suffi- 
cient to heat such small cells up to ~10^ K. 

Second, utilizing the high dynamic range and the Eulerian 
nature of AMR, we have followed the evolution of the hot 



diffuse medium of gaseous halos and the IGM as far as 1 
Mpc away from the galaxies. This allows us to explore the 
interplay between the galactic outflows and the embedding 
medium and to demonstrate that a hot meta l-rich ha lo forms 
around the galaxies from stellar feedback (Cox et al. 2006a). 
The massive gas expulsion in low mass merging galaxies 
leads to a high mass-to-light ratio, as it creates a merger rem- 
nant without much cold gas left for later star formation. 

Although it should be considered provisional, our result 
brings compelling evidence that AMR delivers a uniquely 
powerful tool in understanding merging galaxies, while it ad- 
dresses several issues SPH has suffered from. Comprehensive 
parameter studies should follow, especially in the efficiency 
of stellar feedback and the metal yields of stars; the results 
should be compared and calibrated with observations such as 
the mass- metallicity rel ation ( Trem on ti et a l. 2( )04), galactic 
outflows (|Martinl |2006|) . galactic mor phology (IPark & Choil 
,2008), and gas to stellar mass ratio (iGavazzi et al.l |2008|) . 
Physics such as UV photoelectric heating, cosmic rays, and 
magnetic fields are missing in this work, but will need to be 
considered in the future. 
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